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ABSTRACT 

The correlation between source galaxies and lensing potentials causes a systematic 
effect on measurements of cosmic shear statistics, known as the source-lens clustering 
(SLC) effect. The SLC effect on the skewness of lensing convergence, S3, is examined 
using a nonlinear semi-analytic approach and is checked against numerical simula- 
tions. The semi-analytic calculations have been performed in a wide variety of generic 
models for the redshift distribution of source galaxies and power-law models for the 
bias parameter between the galaxy and dark matter distributions. The semi-analytic 
predictions are tested successfully against numerical simulations. We find the relative 
amplitude of the SLC effect on 53 to be of the order of five to forty per cent. It depends 
significantly on the redshift distribution of sources and on the way the bias parameter 
evolves. We discuss possible measurement strategies to minimize the SLC effects. 
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1 INTRODUCTION 



Recent detections of the cosmic shear signal have opened a 
new window to probe the distribution of matter in the Uni- 
verse, its evolution, and to test cosmological models (Van 
Waerbeke et al. 2000; Wittman et al. 2000; Bacon, Refregier 
& Ellis 2000; Kaiser, Wilson & Luppino 2000; Maoli et al 
2001). These detections have been obtained from relatively 
small fields so far, which limits the statistical analysis of the 
surveys to second order moments, the variance or two-point 
correlation function of cosmic shear. The amplitude of sec- 
ond order statistics reflects that of density fluctuations and 
roughly scales as a 1 oc f2^ 6 ~°' 8 <78 at large scale (Bernardeau, 
Van Waerbeke & Mellier 1997, hereafter BvWM97) and 
cr 7 oc f^ 6 " ' 8 ^' 2 " 1 ' 3 at small scale (Jain & Seljak 1997; 
Maoli et al. 2001). On the other hand, the skewness (a third 
order statistic) of lensing convergence is known to be sensi- 
tive to f2 m , almost independently of <rg (BvWM97). There- 
fore, combined analysis of the skewness and the variance will 
provide precious constraints on both values of fi m and 0%. 
As a consequence, skewness detection and measurement is 



one of main goals of on-going wide field cosmic shear surveys 
such as the DESCART project^. 

Cosmic shear statistics have been studied analytically 
(see Mellier 1999 and Bartelmann & Schneider 2001 for re- 
views and references therein) as well as numerically (Jain, 
Seljak & White 2000; White & Hu 2000). The skewness of 
lensing convergence was first calculated by BvWM97 based 
on a quasi-linear perturbation theory approach. It has been, 
however, recognized that this approach is not robust enough 
to provide accurate predictions for the value of the skewness 
over the whole available dynamic range. In particular, the 
two following points have to be addressed and carefully in- 
cluded in the calculations: (i) Nonlinear growth of the den- 
sity field: numerical studies show that nonlinear growth en- 
hances skewness especially at angular scales smaller than 
one degree (Jain et al. 2000; White & Hu 2000; Van Waer- 
beke et al. 2001b). (ii) Source clustering: Bernardeau (1998) 
(hereafter B98) pointed out that correlations between source 
galaxies and lensing potential reduce skewness amplitude. 



For more information a bout DESCART project, see 



http:/ /terapix. iap.fr/Descart/ 
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B98 underlined that this effect is sensitive to the redshift 
distribution of sources. 

The purpose of this paper is to examine the effect of 
source clustering (SLC) on measurements of the skewness of 
lensing convergence. Special attention is payed to its depen- 
dence on the redshift distribution of sources and on evolution 
of the bias relation between matter and galaxy distribution. 
Since the redshift distribution of faint galaxies is uncertain 
and little is known about the bias, this paper does not aim at 
making accurate predictions for the amplitude of SLC effect 
in real cosmic shear surveys. Our objective is to estimate its 
magnitude in order to propose strategies that minimize its 
effects. 

We basically follow the perturbation theory approach 
first developed by B98 but generalize it in two ways: (i) we 
take into account the effects of nonlinear evolution of the 
density field, adopting the nonlinear semi-analytic ansatz 
developed by Jain & Seljak (1997) and Van Waerbeke et 
al. (2001b); (ii) we allow a possible redshift dependence of 
the bias parameter, b(z) = bo(l + z) 1 and examine the cases 
7 = to 2. Moreover we consider three cosmological Cold 
Dark Matter family models (CDM), two flat models with 
and without cosmological constant and an open model, and 
12 different models for the source distribution which cover a 
wide range of mean redshift and width for the distribution. 

Finally, for the first time the accuracy of semi-analytic 
predictions for the SLC effects on the skewness is tested 
against numerical simulations in standard CDM model. 

The outline of this paper is as follows. In section 2, the 
physical mechanism of SLC is described. In section 3, an ex- 
pression for the skewness of lensing convergence is presented 
that takes both SLC and nonlinear evolution of the density 
field into account. In section 4, our models are described. Re- 
sults of the semi-analytic approach are presented in section 
5. In section 6, semi-analytic predictions are tested against 
numerical simulations. We summarize and discuss our con- 
clusions in section 7. The derivation of the convergence skew- 
ness in the presence of the SLC is presented in Appendix A. 
The JV-body data sets and the ray-tracing method used for 
this work are described in Appendix B. In Appendix C, de- 
tails of the procedure to generate mock galaxy catalogues 
are presented. 



2 WHAT IS THE SLC EFFECT? 

The SLC effect discussed in this paper comes to light be- 
cause of the conjunction of three circumstances, namely: 
(i) source galaxies are not randomly distributed in the sky 
but are correlated; (ii) the source galaxy distribution traces 
somehow the matter field; (iii) the redshift distribution of 
source galaxies is rather broad. The width of the distribu- 
tion depends on source selection criterion, and generally, the 
distribution of source galaxies overlaps with the distribution 
of lensing structures, so source galaxies are somehow cor- 
related with the lensing potential. This correlation causes 
systematic effects on measurements of cosmic shear, that 
may be illustrated as follows. Figure [l] shows a distribution 
of sources (denoted by filled circles) and the gravitational 
potential (contour lines). For line-of-sight 1 (LOS 1), the 
distant galaxies are lensed by the gravitational potential lo- 
cated at an intermediate distance and thus have a high posi- 



LOS 2 




Obs 



Figure 1. An illustration of the correlation between the gravi- 
tational potential (contour lines) and the population of sources 
(denoted by filled circles). 



tive lensing convergence^ This high signal is reduced by the 
excess of foreground sources bound to the foreground grav- 
itational potential which, in contrast, has a low lensing. On 
the other hand, for a line-of-sight 2 (LOS 2), distant sources 
are lensed by the foreground void and thus have a nega- 
tive lensing convergence. This negative signal is amplified 
because of the lack of foreground sources in the void. Ac- 
cordingly, the probability distribution function of the lens- 
ing convergence, which is skewed toward a high value in ab- 
sence of the SLC (e.g., BvWM97, Jain et al. 2000) becomes 
more symmetric than for the case of a random distribution 
of source galaxies. As a result, the amplitude of skewness of 
lensing convergence drops. 

As was pointed out by B98, there is another possible 
effect caused by intrinsic clustering of source galaxies. The 
average distance of sources and their sky density may in- 
deed vary from one direction to another, which can cause 
additional systematic effects on the cosmic shear statistics 
(these effects was included in the definition of source clus- 
tering used by B98). It was pointed out by B98 and Thion 
et al. (2001) that this effect on the convergence skewness is 
in general very small. Therefore, in this analysis, we do not 
take it into account for the analytical calculations presented 



t The lensing convergence is not a direct observable but is ob- 
tained via a convergence reconstruction technique (Van Warbeke, 
Bernardeau & Mellier 1999) or the aperture mass (Schneider et 
al. 1998) from a lensing shear map. 
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in next section, although it will be obviously present in the 
numerical experiments discussed in § 6. 



3 THE PERTURBATION THEORY 
APPROACH 

3.1 The quasi-linear regime 

The expressions for the skewness of lensing convergence 
and the correlation term due to SLC were first derived by 
BvWM97 and B98 in the framework of perturbation theory. 
In this subsection, we only summarize expressions which are 
directly relevant to this paper. Detailed derivations and no- 
tations are given in Appendix A. 

In the presence of the SLC, the skewness parameter, de- 
fined by 5*3(0) = (Kg) /Vk(6), consists of two terms; namely, 
one arises from the quasi-linear theory and the other from 
SLC: 



with 



(1) 



sf(0) = 



3\q.l. 



6 fH^ 3 ' v " 



V*{9) 

x 



(?)■ 



(2) 



sf c (e) 



x / dX P(x',xMx') J o(x') 



G_Ho_ 

Me) c j 

In these equations, 



(3) 



• 9 is the size of the smoothing window, 

• Ho/c is the Hubble constant in speed of light units, 

• Vk (9) is the variance of the lensing convergence [given 
by eq. ju§)], 

• x denotes the radial comoving distance (and xh corre- 
sponds to the horizon), 

• w(x) is the lensing efficiency function [eq. (A.6)], 

• functions Iq(x,9) and h(x,9) depend on the p ower- 



spect rum of matter density fluctuations through eqs. (A13) 
and ( |A16| ), 

• fi m is the matter density in units of the critical density, 

• n s [x(z)] the average number dens ity o f sources, 

• function gixuXs) is given by eq. (A5), 

• b(x) is the linear biasing function between the galaxy 
and the matter density contrast. 

3.2 Nonlinear regime 

For the variance of the lensing convergence, the effect of 
nonlinear evolution of the density power spectrum can be 



included by replacing the linear power spectrum (which en- 
ters the above expressions through V K (9), io(x) an d ^i(x)j 
see Appendix A for their explicit expressions) with the non- 
linear power spectrum, i.e., Pn n (a,k) — > Pt^i J (a,k) (Jain & 
Seljak 1997). We use the fitting formula of nonlinear power 
spectrum given by Peacock and Dodds (1996). This semi- 
analytic approach has been tested against ray-tracing simu- 
lations, and a good agreement between the numerical results 
and the semi-analytic predictions was found (Jain et al 2000; 
White and Hu 2000). 

In the framework of perturbation theory, all density 
contrasts needed for the calculation of the skewness correc- 
tion term [equation (|^)] correspond to linear order (see B98 
for details) . This comes from the fact that the quantity to be 
computed, an angular average of projected fourth moments, 
is given by the product of two two-point correlations: its in- 
trinsic connected part has a negligible contribution because 
of the projection effects. The incorporation of the nonlin- 
ear effects is then straightforward. As for the variance, it 
amounts to formally replacing the linear power spectrum 
with the nonlinear one. 

The semi-analytic calculation of the skewness in the 
nonlinear regime was developed by Van Waerbeke et 
al. (2001b). It is based on the fitting formula of the den- 
sity bispectrum by Scoccimarro & Couchman (2000) and is 
given by 



Sf l - (9) = 



V2{9) (2 



6_ AHbV 
tt) 4 V c J 



dx™ (X) 



x / d\ 1 P N i J (k 1 )W[k 1 f(x)9] 



x / d 2 k 2 P NL (k 2 )W[k 2 f(x)9] 



xiy[|k 1+ k 2 |/( X )0]Fr(kx,k2) 



(4) 



with 



F 2 cff (ki,k 2 ) 



7 



a(n, ki)a(n, fe) 



^ i( ,, : ,, f; ,, 2) ^(| + J) 
(ki ■ k 2 ) 2 



+-c(n, ki)c(n,k 2 )- lu l2]k |2 



[fci] 



(5) 



Here, the notations are as follows: 



• k= |k|, 

• fix) denotes the comoving angular diameter distance 
(see Appendix A) , 

• W is the Fourier transform of the smoothing window, 

• functions a(n,k), b(n,k) and c(n, k) depend on the ef- 
fective power spectral index n at scale k (explicit expres- 
sions are given in Scoccimarro & Couchman 2000; see also 
Van Waerbeke et al. 2001b). 

It should be noted that, since the bispectrum fitting for- 
mula is constructed via the dark matter bispectrum mea- 
sured from only one A r -body simulation data set, there is 
about a 10-20 per cent uncertainty in the fitting formula. 
This is mainly a cosmic variance effect (Van Waerbeke et 
al. 2001b). 
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Table 1. Cosmological parameters for the models considered in 
this paper. f2 m is the matter density in units of the critical density 
and similarly for the cosmological constant Q\ ; h = Ho / 100 is the 
Hubble constant in units of 100 km s" 1 Mpe - , ag is the variance 
of the dark matter distribution in a sphere of radius 8 Mpc. 



Model 


S 2 m 


n x 


h 


0"8 


SCDM 


1.0 


0.0 


0.5 


0.6 


OCDM 


0.3 


0.0 


0.7 


0.85 


ACDM 


0.3 


0.7 


0.7 


0.9 



Table 2. Parameters in n s (z). 



Model (z) 



Al 


1.2 


0.572 


2 


1.5 


0.798 


A2 


1.2 


0.456 


3 


1.8 


0.813 


A3 


1.2 


0.297 


5 


3.0 


1.01 


A4 


1.2 


0.182 


8 


6.0 


1.18 


Bl 


1.5 


0.866 


2 


1.0 


0.500 


132 


1.5 


0.618 


3 


1.5 


0.812 


B3 


1.5 


0.400 


5 


2.5 


1.11 


B4 


1.5 


0.244 


7 


6.0 


1.51 


CI 


0.9 


0.429 


2 


1.5 


0.598 


C2 


0.9 


0.342 


3 


1.8 


0.610 


C3 


0.9 


0.240 


5 


2.5 


0.667 


C4 


0.9 


0.136 


8 


6.0 


0.884 



4 MODELS 

4.1 Cold dark matter models (CDM) 

We discuss three Cold Dark Matter (CDM) models, a fiat 
model with (ACDM) and without cosmological constant 
(SCDM) and an open model (OCDM), using galaxy clus- 
ter abundances to normalize the power-spectrum (Eke, Cole 
& Frenk, 1996; Kitayama & Suto 1997) and the formula 
of Bond & Efstathiou (1984) for the transfer function. The 
parameters of the models are listed in Table [l]. 

4.2 Redshift distribution of source galaxies 

We assume that n s (z) takes the form, 
dz/3 



n s {z)dz 



*»r[(i + a)//3] 



(£) 







) exp 


\z* ) 







(6) 



where T(x) is the Gamma function. 

We explore 12 models for the shape of the distribution. 
The parameters in each model are listed in Table |^. The av- 
erage redshift is (z) = 1.2, 1.5 and 0.9 for models Al-4, Bl-4 
and Cl-4, respectively. We characterize the width of the dis- 
tribution by the root-mean-square, Az, which varies within 
a factor of ~3.2, 3.5 and 3.2 in models Al-4, Bl-4 and Cl-4, 
respectively. Note that only model Al matches roughly the 
observed redshift distribution of galaxies in current cosmic 
shear detections (Van Waerbeke et al. 2000). However, to 
keep our approach as general as possible, we still use a rea- 
sonably large parameter range for the possible shapes of the 
distributions. 

Figure |i| shows the redshift distribution of sources and 
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Figure 2. The redshift distributions of sources, n 8 (z), and the 
corresponding lensing efficiency functions divided by the density 
parameter, ui(z)/f2 m , as functions of redshift. Top panel: for three 
source distribution models in SCDM model. Bottom panel: for Al 
model in three cosmologies. 



the corresponding lensing efficiency as functions of redshift. 
In the top panel, SCDM is supposed, and various models for 
galaxy number counts are taken. In the bottom panel, model 
Al is assumed for number counts, and various cosmologies 
are considered. Roughly speaking, the amplitude of SLC is 
controlled by the amplitude of overlapping between the pop- 
ulation of sources [n s (z)] and that of lenses [which is very 
closely related to w(z)/fl m ]. It is important to keep in mind 
that the normalized efficiency function w(z)/f2 m increases 
in order of SCDM, OCDM and ACDM. 



4.3 Model for the bias 

We assume that the bias between the galaxy and the mat- 
ter distribution is linear and takes a power-law form as a 
function of redshift, i.e., 



b( z ) = b (i + z y 



(7) 



We examine three cases, y = 0, 1 and 2, and we shall take 
feo = 1. Since, so far, little is known about a realistic descrip- 
tion of the bias, we adopted this model for its simplicity and 
the wide possible range of possibilities it nevertheless covers. 
Numerical studies of dark matter clustering combined with 
measurements of two-point correlation function in galaxy 
catalogues suggested that bo is close to unity (e.g., Jenkins 
et al. 1998). 
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Figure 3. Upper panel: predicted skewness of the lensing con- 
vergence for Al and SCDM, with and without source clustering 
effect taken into account. Lower panel: the ratio R as defined in 
equation (H). 



5 RESULTS 

Let us introduce the parameter which characterizes the am- 
plitude of the SLC effects defined by 



R = 



^3 

S 3 



(8) 



where Sf c is the correction brought by SLC and 53 is the 
skewness in the case of the absence of SLC (S3' ' and S3' 1 ' 
for quasi-linear and nonlinear computation, respectively). 

Figure ^ shows S3 (upper panel), with and without tak- 
ing the SLC effect into account, and R (lower panel) as a 
function of 9 for the Al model and SCDM. Nonlinear effects 
on the skewness are discussed in detail in Van Waerbeke et 
al. (2001b). It should be noted that nonlinear growth of the 
density field enhances the skewness significantly at scales be- 
low 1 degree, so the SLC correction term remains relatively 
small because of cancellations between the numerator and 
the denominator in the first line of equation ^ (see also 
Appendix A for the explicit expressions for V K and S§ lc ). 
As a consequence, R decreases significantly when 9 < 20-30 
arcmin. It should be also noted that for 9 J; 100 arcmin, 
where nonlinear effects can be safely neglected, SLC effect 
is reduced while 9 increases. This is due to the change in 
the slope of the density power spectrum occurring when the 
spatial smoothing scale f(x)6 at z — 0.3 ~ 0.8, where the 
most of the lensing contribution comes from, is of the order 
of 10(fi m ft 2 ) -1 Mpc. 

Let us now discuss the theoretical predictions that take 
into account nonlinear effects. Figure ^| shows S3 and R for 
three cosmologies (top panel) and three bias evolution mod- 
els (bottom panel). The top panel clearly shows that it is 
essential to take SLC into account to put constraints on val- 
ues of f2 m determined from S3 . Figure ^ also suggests that 
SLC is more important for low than for high density models. 
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Figure 4. 53 and R [equation (ph] with and without the source 
clustering effect taken into account as functions of scale. The non- 
linear ansatz is used. Cosmology and source distribution models 
are denoted in each plot: (a) three different cosmological mod- 
els are considered; (b) three different bias evolution models are 
considered. 



This is explained by the fact that the efficiency function is 
larger in the first than in the second illustrated by 

Figure ^| The bottom panel of Figure ^ also shows that the 
SLC effect increases with strength of evolution in bias with 
redshift. This is a natural consequence of the fact that the 
SLC effect is caused by the correlation between the lensing 
potential (the matter distribution) and the distribution of 
source galaxies. Finally, note that for 9 10 arcmins, the 
relative SLC effect is nearly independent of scale. 

Figures |^ and ^ show that the relative SLC effect, R, 
peaks around 9 =30-60 arcmin. One might wonder what 



© .... RAS, MNRAS 000, 



6 T. Hamana et al. 



o.i = 



0.1 



- , 1 1 1 , 1 1 , 1 1 1 1 1 1 1 - 

■ SCDM 

O ? 

: x A E 
: o a : 

" a A 7=0 

- x * 7=1 

: a O 7=2 : 

- 1 1 1 1 1 1 1 1 - 


1 < < < < < 
■ SCDM 

-So=10' ox 

: x * • 

o a : 

• O a 

f from left to right 
: A Model A4, A3, A2, Al : 
III — 1 — 1 — 1 — — 1 — 1 — 1 — — I— 


SCDM, 7=1 

_e =r o a 

1 O A = 
O 

A x 

a Model Al-4 
_°a x x Model Bl-4 _ 
E o Model CI -4 \ 


SCDM, 7=1 
_6„=10' Oa *_ 

; Oa x i 

° a x ! 

O a x 

from left to right : 
Model 4, 3, 2, 1 I 


- 1 1 1 1 | 1 1 1 | 1 1 1 | 1 '■ 
"7=1 o 

^B =r o x 

: o 2 : 

A 

" ° a SCDM 
_ a x OCDM 
O ACDM 

- , 1 , , , 1 , , , 1 , , , 1 , - 


- 1 | 1 1 1 1 1 1 1 1 1 1 - 
"7=1 g 

_e„=io' 9 a 

\ 9 A 

: ? * 

A 

from left to right 
: Model A4, A3, A2, Al 



0.2 0.4 0.6 0.B 
Az 



0.4 0.6 
Az 



Figure 5. R as a function of Az as predicted from our semi- 
analytic model. The smoothing scales are 1 arcmin and 10 arcmin 
for left and right panels, respectively. Top panels are for three 
bias models in SCDM, Al-4, middle panels are for all 12 source 
distribution models in SCDM 7 = 1 case, and bottom panels are 
for three cosmological models in 7 = 1 Al-4 cases. 



should be the ideal smoothing scale for measuring the skew- 
ness while reducing as much as possible the SLC effect: 
should it be larger or smaller than the peak position? To an- 
swer this question properly, one also has to consider signal- 
to-noise ratio, S/N. Typically, signal-to-noise in 5*3 is ex- 
pected to decrease with 8 due to the finite area covered by 
the survey. Van Waerbeke, Bernardeau & Mellier (1999) nu- 
merically investigated the efficiency of weak lensing surveys, 
taking into account both this effect and the noise due to 
intrinsic ellipticity of the source galaxies. Figures 8 and 10 
of their paper indicate that it might be difficult to detect 
the skewness with S/N > 1 at smoothing scales larger then 
60 arcmin, even with a wide field survey covering 10° x 10°. 
This suggests that the best choice for the smoothing scale, 
keeping both SLC effects low and a good signal-to-noise ra- 
tio, should be 8 of order of 1 arcmin.jj] 

Figure |E] shows R for 8 = 1 arcmin (left panels) and 
10 arcmin (right panels) as a function of the source redshift 
distribution width, Az. A comparison of the left and right 
panels confirms a visual inspection of Fig. ^, namely that R 
is fairly insensitive to 9 in the scaling regime considered, 8 ^ 
10 arcmin. The top and bottom panels indicate that effects 
of cosmology and bias evolution model on the amplitude of 
parameter R are significant, but the shape of R as a function 
of Az remains fairly stable. Note furthermore that in the 
middle panels, models with the same mean source redshift 
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Figure 6. R computed by the semi-analytic formula versus that 
derived from the phenomenological law [eq. (^j)]. The parameters 
(A, B, C) for each case are summarized in Table ^. 



Table 3. Parameters of the phenomenological law [eq. (|9|)] de- 
rived from semi-analytic calculations 
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■f It should be however noted that to break the degeneracy be- 
tween cosmological parameters, one still has to measure cosmic 
shear statistics at linear scales, i.e. 8 > 1 degree (Jain &: Seljak 
1997). 



form sequences in the R-Az plane with very similar slopes 
at fixed Az. 

This suggests that for a choice of the cosmological model 
and 7, there exists a simple phenomenological law that re- 
lates R to (2) and Az, which is valid for all source models 
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<z> 



Figure 7. Contour lines of R obtained from our semi-analytic 
modeling in {z)-Az plane for 9 = 1' and 7 = 1. Each panel 
corresponds to a different cosmology. 



considered here, 
R = A 



0) 



where A is of the order of 0.5, B and C varying from 1.5 to 
3. The precise values of parameters (^4, B, C) are obtained 
by a least-squares fitting method. They are given in Table 

S a 

The accuracy of this law is demonstrated in Figure |6| 
One can see that the data lie fairly well on the parameterized 
line. It allows us to make a contour plot of parameter R 
in {z)-Az space, as shown in Fig. [?], for three cosmological 
models. Not surprisingly, this figure clearly indicates that 
the way to reduce out SLC effect on a measurement of the 
skewness is to make the source distribution narrow with a 
high mean redshift. 

Finally, we examine the dependence of the SLC effect 
on the evolution of bias. Figure M shows R as a function of 
7 for 4 models selected arbitrarily. For each source distribu- 
tion model, the R-y relation is well fitted by the following 
empirical law, 



log R — £7 + const. 



(10) 



The coefficient e, which describes the strength of dependence 
of the SLC effect on the evolution of bias, is shown in Fig. ^ 
as a function of R(y = 1) for all source distribution models 
we consider. Each panel corresponds to a given choice of 
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Figure 8. Semi-analytic values of R as a function of 7 for 4 
selected models in SCDM. 
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model. Coefficient e as defined in equation ( |lo[ ) is shown as a 
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cosmology. One can see that e remains in the range 0.1 < 
e < 0.3 and is almost insensitive to both smoothing scale 
and cosmology. At a fixed value of R(j = 1), e increases 
with the mean redshift (in order of C, A and B). This is a 
natural consequence of the fact that, for our choice of the 
bias evolution model eq. (Q) , the impact of the change in the 
bias evolution is more significant at a higher redshift. 

The uncertainty in parameter R caused by our igno- 
rance of b(z) can be roughly estimated using the empirical 
relation (Flfj) as follows: suppose that the power-law model 
([?]) for the evolution of bias stands, but that there is an er- 
ror A 7 on the value of 7. Applying simple error propagation 
technique, one finds SR/R — 2.3eA 7 (~ 0.5A 7 ): if one is 
able to constrain the bias evolution model with an accuracy 
better than A 7 < 0.4, the uncertainty in R drops below 20 
per cent. 



6 TESTING SEMI-ANALYTIC PREDICTIONS 
AGAINST NUMERICAL SIMULATIONS 

In this section, we compare theoretical predictions to ray- 
tracing experiments in JV-body simulations, using mock 
galaxy catalogues extracted from the simulations as dis- 
tributions of sources. The iV-body data sets and the ray- 
tracing method used for this work are described in Ap- 
pendix B. In Appendix C, details of the procedure to gen- 
erate mock galaxy catalogues are presented. We focus only 
on the SCDM model, but the conclusions of our numerical 
analysis should not depend significantly on the considered 
cosmology. 

We measure convergence statistics on mock galaxy cat- 
alogues as follows: the value of convergence for a galaxy at 
redshift z is given by linear interpolation between the m 
— computed from rays propagating between redshift z of 
the closest lens plane to the galaxy (see Appendix B) and 
present time — measured at the four nearest angular pixels 
from the galaxy position in the sky. The amplitude of the 
SLC effect is measured by comparing simulations to similar 
SLC-free ones. They are obtained from other mock galaxy 
catalogues with the same source distribution n 3 (z) (i.e. the 
same as in the SLC mock catalogue), in which galaxies are 
randomly distributed on the sky. Note finally that top hat 
filtering is used and the intrinsic ellipticity of galaxies is not 
taken into account. 

The upper panel of Fig. [j^ shows the skewness pa- 
rameter measured from the 7 = mock galaxy catalogue 
with and without the SLC effect. The lower panel displays 
the function R(0), except that in the denominator of equa- 
tion (B), we always take the value given by nonlinear semi- 
analytic predictions, S3". As discussed in Appendix B (see 
also Hamana & Mellier 2001), the simulations have lim- 
ited available dynamic range since they are contaminated 
by force softening and finite volume effects at small and 
large scales respectively, where the measured S3 is expected 
to underestimate the real value. Furthermore, there is a 10- 
20 per cent uncertainty in the nonlinear perturbation theory 
predictions. With these elements in mind, we see that agree- 
ment between measurements and predictions is reasonable 
when SLC effects are taken into account. In particular, the 
order of magnitude of the shift between the upper and the 
lower symbols in the top panel of Fig. [u^ matches very well 
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Figure 10. 53 (upper panel) and R (lower panel) as functions of 
the smoothing angle 9 in the 7 = case. The solid (dotted) line 
shows the semi-analytic prediction with (without) the SLC effect 
taking into account, while the filled circles (triangles) show the 
results of ray-tracing simulations with (without) SLC effect taking 
into account. Error bars give some estimate of the uncertainty E 
on the measurements, i.e. E(Si) ~ AS3 / v / 40, where (AS3) 2 is 
the dispersion over the 40 realizations. Note that R is calculated 
as -R = — Sg /S3 ' 1 where S3 ' 1 ' is always given by semi-analytic 
predictions. 

that between the dotted and the solid curve, as illustrated 
by bottom panel. 

Similar results are obtained for the 7 = 1 and 2 cata- 
logues, as summarized in Fig. jH], which concentrates on the 
parameter R. Since numerical experiments for 7 = 1 and 2 
cases are done with threshold bias instead of the linear bias 
used in the semi-analytic computation, we should only fo- 
cus on the differences among the three models arising from 
different evolution of biasing. Although there is a system- 
atic difference between the predictions and measurements, 
a trend in the dependence of evolution of biasing on the 
SLC effect found in the predictions is well reproduced in the 
measurements. 

Furthermore, one finds that the 7 = 2 measurement 
agrees better with the 7 = 3 prediction except for the largest 
angular scales. This is consistent with the scale dependence 
of biasing detected in the 7 = 2 mock catalogue (on small 
scales the bias evolves as 7 — 3, see Appendix C for de- 
tailed discussion on this point). We may conclude from the 
results above that the semi-analytic approach gives a good 
prediction of the SLC effect on the convergence skewness. 



7 SUMMARY AND DISCUSSION 

We have examined the source-lens clustering (SLC) effect 
on measurements of the skewness of lensing convergence us- 
ing a nonlinear semi-analytic approach. The result of semi- 
analytic predictions were tested against numerical simula- 
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Figure 11. The parameter R as displayed in the lower panel of 
Fig. but now for all values of 7 considered. Symbols show the 
results of ray-tracing simulations: circles, triangles and squares 
correspond respectively to 7 = 0, 1 and 2 mock catalogues. Lines 
correspond to semi-analytical predictions, 7 = 0, 1, 2, 3 from 
bottom to top. Error bars are computed as explained in caption 
of Fig. |h]. 

tions, and a good agreement between them was found. Our 
main conclusions are as follows: 

• SLC effect strongly depends on the redshift distribution 
of source galaxies. We found that the effect scales with the 
width and mean redshift of the distribution roughly as R oc 
^-(3.0-1.8)^1.4-2^ ,- where R = sf c /S 3 , and Sf c is the 
change in measured S3 due to SLC). As illustrated by Fig. ^ 
this relation indicates that it is essential to make the width of 
the distribution narrow and its mean redshift high to reduce 
the SLC effect (this was partly pointed out by B98). 

• SLC effect also depends on the evolution of the bias 
between the galaxy and total matter distributions, b(z). As- 
suming a simple power-law model and linear bias, b(z) oc 
(1 + z) 7 , we found that the uncertainty in 7 transforms into 
8R/R — 2.3 eA 7 with a typical value of e ~ 0.2. This in- 
dicates that the uncertainty in 7 must be A 7 < 0.2 for 
predicting the amplitude of the SLC effect with better than 
10 per cent accuracy. 

The main uncertainty in semi-analytic predictions 
comes from the fact that the accuracy of the nonlinear fit- 
ting formula of the density bispectrum is only 10-20 per cent. 
We expect the same level of uncertainty in predictions pre- 
sented in this paper for the SLC effect on the convergence 
skewness. This can actually be improved by measurements 
in large A-body simulations with high spatial resolution. 

Since, so far, little is known about the evolution of the 
bias, it is still very difficult to predict the SLC effect accu- 
rately. It is therefore very important to reduce this effect 
as much as possible by controlling the redshift distribution 
of sources. The above results tells us that an ideal obser- 



vational strategy might be as follows: (i) going to a deep 
limiting magnitude to increase the mean redshift of the sur- 
vey and (ii) using only fainter images to reduce the width of 
the distribution. A desirable source distribution for R < 0.1 
suggested by Fig. ^ would have Az < 0.3 and (z) > 1. 
This may be of course challenging: going to deeper mag- 
nitudes will make the calculation of the redshift distribu- 
tion of sources more difficult, and using only faint images 
will increase the noise due to intrinsic ellipticity of galaxies. 
We leave more detailed studies on the designing of optimal 
strategies to future works. 

Before closing this paper, it is important to empha- 
size that the use of source galaxies with a small redshift 
width may unfortunately introduce additional skewness sig- 
nal due to the intrinsic correlation of galaxy ellipticities. It 
was indeed suggested that the amplitude of this latter ef- 
fect scales roughly as (Az)" 1 , although the normalization 
of this relation is ambiguous because of the uncertainty in 
the correlation between the shape of galaxies and that of 
their dark matter halos (Croft & Metzler 2000; Crittenden 
et al. 2000; Heavens, Refregier & Heymans 2000; Pen, Lee 
& Seljak 2000; Catelan, Kamiokowski & Blandford 2001; 
Hatton & Ninin 2001). However, as pointed out by Croft & 
Metzler (2000), intrinsic ellipticity correlations might just 
act as an additional source of random noise, without signifi- 
cantly influencing the measured value of the skewness of the 
convergence. 
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APPENDIX A: PERTURBATION THEORY 
APPROACH TO THE COSMIC SHEAR 
STATISTICS IN THE PRESENCE OF SLC 

The expressions for the skewness of lensing convergence 
and the correlation term due to SLC were first derived by 
BvWM97 and B98, respectively, in the framework of per- 
turbation theory. However, B98 only gave the expression for 
the case of an Einstein-de Sitter cosmological model and as- 
sumed power-law density power spectrum. In this Appendix, 
for this paper being self-contained, we re-derive the skewness 
terms which are valid for an arbitrary Friedmann model. 
We basically follow the original deviation by B98. It should 
be noted that the skewness corrections are not only caused 
by the SLC but also arise from e.g., the lens- lens coupling 
(BvWM97, Van Waerbeke et al 2001b) and the lensing mag- 
nification effect (Hamana 2001). In what follows, we focus 
on the SLC correction and are not concerned with other 
correction terms. 



Al Fluctuation in a source distribution due to 
the source clustering 

The number density of the sources at redshift z and in di- 
rection d> can be written, 



< bs [xM,</>] = n s [ X {z)] (1 + S a \x(z), 0]) , 



(Al) 



where n 3 [x(z)] is the average number density of sources, 
5 s [x(z),4>] is their local density contrast and x denotes 



the radial comoving distance. We suppose as B98 that 
the average source number density is normalized to unity, 
f* H dxn s (x) ~ lj where %h is the distance to the horizon 
(the normalized distribution denotes the probability distri- 
bution). Following B98, we assume that the density contrast 
of sources is related to the matter density contrast, 5, via 
the linear biasing, 



S s (x,<P) = b(x)5(x,<P)- 



(A2) 



A2 Convergence statistics in the presence of SLC 

Let us consider the measured convergence that results from 
averages made over many distant galaxies located at differ- 
ent distances. Denoting the smoothing scale by 9, such an 
average can formally be written as, 

E£i*W>.(*.*i) (A:!) 



where Wg(x) denotes the weight function of the average, N„ 
is the number of source galaxies, K 3 (zi, 4 ) is the lensing con- 
vergence signal from a source galaxy located at redshift Zi in 
a direction i and is given by (e.g., Mellier 1999; Bartelmenn 
& Schneider 2001 for reviews) 



K a {z,4>) = 



with 



3f2 m Ho 



%.{*) 



dXig{XhXs)8{Xh^), 



g(x.hX.s) = 



Ho f(xi)f(Xs - Xi) 



(A4) 



(A5) 



c f(Xs)a(xi) 
Here a is the scale factor normalized to its present 
value, and f(x) denotes the comoving angular diame- 
ter distance, defined as fix) = K -1 ^ 2 sin K 1 ^ 2 x, X> 
(-K)- 1/2 sinh(-K) 1/2 x for K > 0, K = 0, K < 0, re- 
spectively, where K is the curvature which can be expressed 
as K = (Ho /c) 2 (fi m — 1). For the weight function, the 
angular top-hat filter (BvWM97) and/or compensated fil- 
ter (Schneider et al. 1998) are commonly adopted (e.g., Van 
Waerbeke et al. 2001a). In what follows, we consider the top- 
hat filter for the weight function, and in this case equation 



(A3) is reduced to k$ — J^^f, n s (zi, 4>i)/N^, where iVj is the 



number of source galaxies within an aperture 8 centered on 
a direction <p-. The number density of sources for current 
and future weak lensing analyses is about 40 per arcmin 2 
(e.g., Van Waerbeke et al. 2001a), which typically implies 
more than 100 galaxies in discs of radius 9 > 1 arcmin. As a 
result, discreteness effects from the source distribution can 



be neglected (see also B98) and we can rewrite (A3) in the 
continuous limit: 



J d 2 <t>w e {<t>) J o XH dx^(x,0K bs (x,</>) 

J d^W 6 (4>) J Q XH d X nt B {x,4>) 



(A6) 



Let us now expand equation ( A6) in terms of 5 using the per- 
turbation theory approach, following BvWM97. The pres- 
ence of SLC does not change the expression of the first order 
term, 



J 1 ) 



3Sl m H[ 



dXs n s (Xs 



dXi9(Xi,Xs)S { >(xi,<t>) 
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with 



where w(x) is so-called the lensing efficiency function defined 
by 



w{xi 



(A8) 



The second order convergence consists of two terms: one 
comes from the second order density perturbation and it is 
formally written by r eplacing the subscript in the first 
order expression @ with (2) (BvWM97); the other one is 
due to SLC, 



sic(2) 3f2 m Ho i j 
k„ w = — ; / d 4>We(4>) 



2 c 

x dxn s (x)b(x)S W (x^) 
Jo 



x / dx'g(x',x)S (1) (x',& 
Jo 

-kP j d 2 4>We((j}) 

x [ " d X n s (x)Kx)5 m (x,<t>)- 



(A9) 



Using the small angle approximation (Kaiser 1992), equation 
(A.7) is rewritten in terms of the Fourier transform of the 



density contrast, 5(k), as 



Ho 



j£k_ 

(2tt) 3 



xW[f( X )k±0], 



5 (1) [fc; X ]exp[ifc x /( X )] 
(A10) 



where the wave vector k is decomposed into the line-of-sight 
component k x and its perpendicular, fej_, and W(x) is the 
Fourier transform of the weight function. In the case of the 
top-hat filter, W(x) = 2Ji(x)/x where Ji is the B esse l func- 
tion of first order. In the same manner, equation ( A9) reads 



K sic(2) _ 30n-Ho j d X n 3 (x)b(x) I dx'g(x',x) 



2 c 



x 



d A k 

d 3 k' 



S w [k;x]eMik x f(x)] 
5 (1) [fe'; X ']exp[ifc' x /( X ')] 



(2tt) 3 

xW[\f( X )k± + f( x ')k\\e] 

rxH 

I dxn s (x)b(x) 
Jo 

d 3 k 



^6 (1) {k;x}eMik x f(x)}W[f(x)ki 



(27T) 



(All) 



The average of the convergence is not affected by the 
presence of the SLC and is therefore zero, («) = 0. The 
variance is not affected by it either at linear order and is 
given by 



V K {6) 



<«r> 



dxi w 2 (.Xi)Io(Xi,Q), 



(A12) 
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Hx,o) 



2tt 



dkkP lin (x,k)Wi [kf(x)6], 



(A13) 



where P iin [x(z),k} is the linear density power spectrum. In 
the presence of SLC, the skewness parameter, defined by 
S3 (9) = (kI) /Vk(8), consists of two terms: one comes from 
the second order perturbation (BvWM97), 



3\ q.l. 



3{fvg ' Kg 



= 6 



dx» (x) 



/_/ , • ! '' x " 

( 

po(x,0) + \lo(x,0)h(x,( 



where 

The other arises from SLC, 



(A14) 



(A15) 



(A16) 



/ 3 \ sic 

(««) 



= 3(4 1)2 4 c(2) > 

= 9n n 



(^°) / d X n s (x)&(xMx)Mx) 



3 /"Xff 



x / d x ' g(x',x)w(x')lo(x 

'0 

tlo 



-6v K (ey 



dxns(x)Kx)w(x)Io(x)- 



(A17) 



To derive the last expression, we used an approximation, 
which turns out to be very accurate for top-hat smoothing 
(see B98 for details), 

f27T 



271 



dtfsiniW(|fc + fc'|) ~ W(k)W(k'), 



(A18) 



where 1? is the angle between the wave vectors k and k' and 
k = \k\. 

Then, the convergence skewness simply reads, in the 
second order perturbation theory framework, Ss(9) = S3' 1 ' + 
Sf c = (Ke> q ' 1 7K 2 (6 , ) + (Ke) sl 7K 2 (6 l ), where all the terms are 
computed above. 

Note that our calculation is slightly different from that 
of B98. Indeed B98 assumed an optimally weighted estima- 
tor for the convergence leading to eq. (9) of his pap er in - 
stead of our eq. (A.6). However, with approxim ation (A18), 

Llll) would 



both estimators give the same results: eq. 
match eq. (22) of B98, and therefore we would easily recover 
eq. (29) of B98 for a scale-free power-spectrum^. 



APPENDIX B: A BRIEF DESCRIPTION OF 
THE RAY-TRACING SIMULATIONS 

In this Appendix, we describe the iV-body data sets and 
the ray-tracing method used for this work. More technical 
details are presented in Hamana & Mellier (2001). 



S Notice the difference of sign convention we use for re, to enforce 
positively for the convergence skewness. 
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compute the skewness of the lensing convergence for single 
source planes, i.e., n s (z) = 5d(z = z s ) where Sd is the Dirac 
delta function. At this stage, we do not take into account 
the SLC effect. Figure 



Bl 



10 15 
0[arcmin] 



Figure Bl. The skewness 53 of the lensing convergence measured 
from the ray-tracing simulations (symbols) compared with non- 
linear predictions (solid lines) for single source redshifts z s ~ 1 
(top), 2, and 3 (bottom). 



Light ray trajectories are followed through large N- 
body simulations data set generated with a fully vectorized 
and parallelized Particle- Mesh (PM) code. Each iV-body ex- 
periment involves 256 2 x 512 particles in a periodic rectan- 
gular box of size (L, L, 2L). The mesh used to compute the 
forces was 256 2 x 512. A light -cone of the particles was ex- 
tracted from each simulation during the run as explained 
in Hamana, Colombi & Suto (2001). Our aim was for the 
light-cone to cover a large redshift range, < z 3, and 
a field of view of 5 x 5 square degrees. To do that, we 
adopted the tiling technique first proposed by White & Hu 
(2000): we performed 11 independent simulations covering 
adjacent redshift intervals [zf 111 , zf lax ], i = 1, . . . , 11. The 
size of each simulation is such that the portion of the light- 
cone in [z™ m , z™ ax ] (aligned with the third axis) exactly fits 
the box-size. This way, angular resolution is approximately 
conserved as a function of redshift, except close to the ob- 
server. Finally, in order to have enough structures in each 
box, we impose the supplementary constraint L > 8O/1 -1 
Mpc. As a result, L follows the following sequence with red- 
shift, 80, 80, 80, 80, 80, 120, 160, 240, 320, 480, 640 h' 1 Mpc. 

The multiple lens-plane algorithm was used for ray- 
tracing calculations (Jain et al. 2000 and references therein). 
The lens planes (which are, at the same time, source planes) 
are separated by intervals of 80h~ Mpc, amounting to a to- 
tal number of 38 in the redshift range < z 3. For each 
ray, the lensing magnification matrix is computed on the 
source planes and is stored. We performed 40 realizations of 
the underlying density field by random shifts of the simu- 
lation boxes in the (x, y) plane. For each realization, 512 2 
rays are traced backward from the observer. The initial ray 
directions are set on 512 2 grids, which correspond to pixels 
of angular size 5°/512 ~ 0.59 arcmin. 

Before using realistic redshift distribution of sources, we 



shows S3 obtained from the sim- 
ulations compared to nonlinear predictions. Measurements 
match theory reasonably well, as expected (Van Waerbeke 
et al. 2001b). There are slight differences which can be ex- 
plained as follows: 

(i) The iV-body simulations have a finite spatial resolu- 
tion, which implies a flattening of 5*3 at scales smaller than 
about 4 arcmin. 

(ii) At large angular scales, 9 ^ 20—40 arcmin, depending 
on the source redshift considered (Hamana & Mellier 2001), 
the measured S3 underestimates the real value, due to finite 
volume effects (i.e. the lack of the large scale power, which 
contributes to the skewness on smaller scales, due to the 
finite size of the simulation boxes, e.g., Colombi, Bouchet & 
Schaeffer 1994; see also Seto 1999). 

(iii) There is an uncertainty in the fitting formula of the 
density bispectrum (section 3), which transforms into a 10- 
20 per cent error on the semi-analytic prediction for S3. The 



differences between theory and measurements in Fig. Bl are 
smaller than this expectation, at least in the range where 
the measurements are reliable, 4 9 J~ 20 — 40 arcmin 
(derived from the above discussion on spatial resolution and 
finite volume effects) . It is important to note this range is not 
equal to the dynamical range that the ray-tracing simulation 
originally have, which is much wider (see Hamana & Mellier 
2001 for discussion on this point). 

In conclusion, without SLC effects (yet) taken into account, 
the semi-analytic prediction obtained for S3 is accurate (Van 
Waerbeke et al. 2001b). 



APPENDIX C: PROCEDURE TO GENERATE 
MOCK GALAXY CATALOGUES 

We generated three mock galaxy catalogues with galaxy 
number counts n s (z) derived from the semi-analytic model 
by Devriendt & Guiderdoni (2000) and reproducing as well 
as possible the power-law model for the function b(z), b(z) = 
(1 + z) 7 with 7 = 0, 1, and 2. They are extracted from ex- 
actly the same dark matter distributions as the ones used 
for the ray-tracing simulations. 

The procedure to create the mock catalogues can be 
described as follows: 

(i) We adopt threshold biasing, i.e. for a smooth density 
distribution of dark matter, we assume that galaxies lie in 
regions of density contrast larger than some threshold which 
may eventually depend on redshift [point (ii) below], 8 > 
Sth(z). Inside these regions, the local number density of 
galaxies at position (z,9,<f)), (where z denotes the redshift 
and (9, (f>) denote the direction in the sky) is proportional to 
dark matter density 



n,(*,M) = A*(*)[l + *(*,M)]- 



(CI) 



The normalization factor fi(z) is such that the redshift dis- 
tribution of galaxies reproduces (in terms of an ensemble av- 
erage) some prior, n s (z), discussed in (iii). To estimate the 
local density contrast from our discrete dark-matter parti- 
cle distribution, we use local adaptive smoothing: the mean 
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quadratic distance d between each simulation particle and 
its 6 nearest neighbors is computed, d 2 — d 2 where di 

is the separation between a central particle and i-the nearest 
neighbors. Then, 1 + 5 oc d~ z . For each dark matter particle 
in regions with 5 > 5th, N galaxies are randomly placed in 
the sphere of radius d centered on the particle position. N 
is computed from a random realization of a Poisson distri- 
bution with average N = (4/3)7rd 3 n g , where n g is the mean 
number density of galaxies. 

(ii) Function 5th (z) is determined numerically so that 
the measured variances of density fluctuations in a sphere 
of radius 8/i _1 Mpc in the galaxy and the dark matter dis- 
tribution, respectively a| aL and as, satisfy 



T gal- 



CO 



= Kz). 



(C2) 



To do that, we use snapshots of the simulations at vari- 
ous redshifts Zi, and compute iteratively 5th (zt) to match 
equation (C2) within a 3 per cent accuracy. Then function 
5th (z) is obtained by linear interpolation of the 5th (zi)- 
Note that since for 7 = 0, b = 1 irrespective of the redshift, 
the galaxy distribution directly traces the matter distribu- 
tion, and thus we do not need to have threshold i.e., we 
simply set 5th = 5 m i n — 1. 

(iii) A prior function n s (z) is needed to compute the nor- 
malization factor fi(z) equation (|Cl[). We have used the 
ab-initio semi-analytic approach to galaxy formation de- 
scribed in Devriendt & Guiderdoni (2000) to obtain a rea- 
sonably realistic estimate of this function. Such an approach 
is based on a Press-Schechter like prescription to compute 
the number of galaxies as a function of redshift, coupled 
to spectro-photometric evolution of stellar populations to 
calculate their luminosities. The results naturally match ob- 
served galaxy number counts and redshift distributions, as 
well as the diffuse extragalactic background light for wave- 
lengths ranging from the UV to the near IR. Here, we sup- 
pose that galaxies are selected in the I band, down to the 
magnitude Iab = 24.5. As a result, the final mock cata- 
logues yields a typical surface number density of 29 sou rces 
per arcmin 2 distributed in redshift as shown in Figure |Cl[ 
The distribution has a peak at z ~ 0.4, with mean redshift 
(z) ~ 0.8 and typical width Az ~ 0.6. 

Note that, for b > 1, the linear bias prescription breaks 
down, because 5 S can not be less than —1 by definition (in 
the perturbation theory approach, this does not cause a se- 
rious problem because |5| is supposed to be much less than 
unity) . Therefore, we do not take the linear bias but use the 
threshold bias for the cases of 7 = 1 and 2. Accordingly, 
strictly speaking, the direct comparison between the semi- 
analytical prediction and the numerical simulation makes 
sense only for the 7 = model. However, we take 7 = 1 
and 2 models to test the effect of the redshift evolution of 
bias which should not depend strongly on details of the bi- 
asing prescription. Moreover, to test the robustness of semi- 
analytic predictions in which the simple linear bias is used, 
it is interesting to use a different biasing prescription for the 
numerical experiments. 

Figure C2 shows the scaling behavior of the bias factor 
defined by 



b(z,£) 



Cgal.(l,*) 

a(£,z) ' 



(C3) 




Figure CI. The distribution of sources (histogram with error 
bars) and the lensing efficiency function (dotted line) as func- 
tions of redshift. Error bars denote standard deviation computed 
among 40 realizations as discussed in Appendix B. 



as measured in the mock catalogues with 7=1 and 7 = 2. 
In this equation, c gal (£) and a (£) are respectively the vari- 
ances in a sphere of radius £ of the galaxy and the matter 
density distribution. In fact, we take for a 2 the variance mea- 
sured in the mock catalogue with 7 = which is unbiased 
by definition. To correct for variations of the selection func- 
tion we use the method proposed by Colombi, Szapudi & 
Szalay (1998). The curves on each panel correspond to red- 
shift slices of [0.16,0.20], [0.3,0.36], [0.56,0.68], [1.09,1.41] 
and [2.40,3.00]. 

By construction, the value of b measured at £ = 8 _1 Mpc 
(triangles) matches very well relation (|C2|). However there 
is no guarantee for this result to hold at all scales. In other 
words, at fixed z, function b(£, z) is not necessarily a constant 
of scale [and equal to b(z) — (1 + z) 7 ], although this is pretty 
much the case for the 7 = 1 mock catalogue. 

For the 7 = 2 mock catalogue, function b(£, z) presents 
large variations with scale, increasing with redshift. This can 
be modeled as a varying effective j(£), for examp le 7 ~ 3 for 
I — 1/t -1 Mpc (squares on bottom panel of Fig. C2). While 
converting scales to angles, more relevant to our analysis, the 
modeling in terms of a function y(0) is not very convincing. 
Still, we find that in the range of interest, 4 ^ 8 ^ 20 — 30 
arcmin, we should compare measured SLC effects to semi- 
analytic predictions corresponding to 2 ^ 7 S~ 3. 



This paper has been produced using the Royal Astronomical 
Society/Blackwell Science IATgX style file. 
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Figure C2. The biasing function b [equation ( JC3| )] as a function 
of spatial scale £, measured in the mock galaxy catalogues. The 
top and bottom panels correspond respectively to 7 = 1 and 7 = 2 
(by construction, we exactly have 6=1 for 7 = 0). Each curve is 
for a fixed value of z, namely z ~ 0.18, 0.33, 0.62, 1.25 and 2.70 
from bottom to top of each panel. On the top panel, the squares 
and the triangle give the values expected from b(z) = 1+z. On the 
bottom panel, the square and the triangle correspond respectively 
to b(z) = (1 + z) 3 and 6(2) = (1 + z) 2 . 
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